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In this work, we consider the possibility that wormhole geometries are sustained by their own 
quantum fluctuations, in the context of modified dispersion relations. More specifically, the energy 
density of the graviton one-loop contribution to a classical energy in a wormhole background is 
considered as a self-consistent source for wormholes. In this semi-classical context, we consider 
specific choices for the Rainbow's functions and find solutions for wormhole geometries in the cis- 
planckian and trans-planckian regimes. In the former regime, the wormhole spacetimes are not 
asymptotically flat and need to be matched to an exterior vacuum solution. In the latter trans- 
planckian regime, we find that the quantum corrections are exponentially suppressed, which provide 
asymptotically flat wormhole geometries with a constant shape function, i.e., b(r) = rt, where r t 
is the wormhole throat. In addition to this analysis, we also fix the geometry by considering the 
behaviour of a specific shape function through a variational approach which imposes a local analysis 
to the problem at the wormhole throat. We further explore the respective parameter range of the 
Rainbow's functions, and find a good agreement with previous work. 

PACS numbers: 04.60.-m, 04.20. Jb 



I. INTRODUCTION 



Modified Dispersion Relations (MDR) are a distortion of the spacetime metric at energies comparable to the Planck 
energy. A general formalism, denoted as deformed or doubly special relativity [l[ was developed in in order to (i) 
preserve the relativity of inertial frames, (ii) maintain the Planck energy invariant and (iii) impose that in the limit 
E/Ep — > the speed of a massless particle tends to a universal constant, c, which is the same for all inertial observers. 
In particular, in curved spacetime the formalism imposes that the relationship between the energy and momentum of 
a massive particle m in special relativity are modified at the Planck scale [HQ], i.e., E 2 g 2 (E / Ep) — p 2 g 2 (E / Ep) = m 2 , 
where the two unknown functions g\(E/Ep) and g%{E / Ep), denoted as the Rainbow's functions, have the following 
property 

lim g 1 (E/E P ) = l and lim g 2 (E/E P ) = l. (1) 

E/Ep^O E/Ep^O 

One interesting aspect of Rainbow's functions is the possibility of regularizing the divergent behaviour of some field 
quantities such as the energy density. This has been done in Ref. [3] with the following choice for the Rainbow's 
functions 

" E i ( E 2 \ 

g 1 (E/Ep) = J2/3 i ^exp^-a^j, g 2 (E/E P ) = l; a>0,ft6i (2) 

More specifically, in Ref. Q, motivated by the promising results obtained in the application of Gravity's Rainbow 
to black hole entropy it was shown that for an appropriate choice of the functions g\{E/EP) and g 2 {E/EP), 
the UV divergences of the Zero Point Energy disappear. It is interesting to note that every choice for the Rainbow's 
functions is restricted by convergence criteria, namely a pure polynomial cannot be used without the reintroduction 
of a regularization and a renormalization process. 
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In fact, MDR have been considered in a wide variety of contexts, such as in cosmology [5] and black hole physics 
|. For instance, the rainbow version of the Schwarzschild line element is given by 

, 2 / 2MG(0)\ dP df 2 r 2 , ... 

ds 2 = - 1 — oTttTT^ n + "7 7T\ + 9/r,/r, \ (dd + sm 2 6»d0 2 ) . 3) 

V r g\{EEp) i_^°b( £ / £p ) &{EE P y 



In particular, note that the position of the horizon, in fixed and energy independent coordinates, is at the usual fixed 
coordinate 2MG(0). However, it was shown that the area of the horizon is energy dependent 0]. 

In this paper, we use the scheme given by Rainbow's function([2"j). in order to verify if the MDR provide interesting 
solutions for wormhole geometries. The latter are hypothetical tunnels in spacetime @, S, an d have been studied in 
a plethora of environments such as dark energy models |j| , modified theories of gravity , observational signatures 
using thin accretion disks [TTj] , and in the semi-classical regime [l2T - [l4j |. amongst many other contexts. In the semi- 
classical regime, the possibility that these wormholes are sustained by their own quantum fluctuations was considered 
in [l3l flij . The graviton one loop contribution to a classical energy in a wormhole background was taken into 
account, and a variational approach with Gaussian trial wave functionals was used. A zeta function regularization 
was involved to handle the divergences and a renormalization procedure was introduced. Thus, the finite one loop 
energy was considered as a self-consistent source for traversable wormholes. The aim of this work is to consider 
the possibility that wormhole geometries be sustained by their own quantum fluctuations, in the context of modified 
dispersion relations. 

This paper is organized in the following manner: In Section [II] we briefly outline the formalism of the classical term 
in Rainbow's Gravity and the one loop energy in a spherically symmetric background. In Section IIII1 we consider 
specific Rainbow's functions and find solutions of self-sustained wormhole geometries in the context of MDR. In 
Section ITVl we fix the geometry by considering the behaviour of a specific metric function b(r), through a variational 
approach, at the wormhole throat, r t . In Section [Vj we conclude. 



II. THE CLASSICAL TERM IN RAINBOW'S GRAVITY AND THE ONE LOOP ENERGY IN A 

SPHERICALLY SYMMETRIC BACKGROUND 

Consider a wormhole spacetime described by the following line element 

ds 2 = -N 2 (r) ? dt ] - + ^ + - lr -^ (dd 2 + sin 2 9dtf) , (4) 

V ' 9l (E/E P ) ^ _ g 2 {E/Ep) g 2 (E/Ep) [ V > > 1 ' 

where N(r) is the lapse function and b(r) is denoted the shape function, as can be shown by embedding diagrams, it 
determines the shape of the wormhole [7]. For simplicity, we consider N(r) = 1 throughout this work. A fundamental 
property of traversable wormholes is the flaring out condition of the throat, given by (b - b'r)/b 2 > [7]. Note that 
the condition 1 — b/r > is also imposed. To be a wormhole solution the following conditions need to be satisfied at 
the throat: b(r t ) — r t and b'(r t ) < 1; the latter follows from the flaring out condition. Asymptotic flatness imposes 
b(r)/r — > as r — > +00. However, one may also construct solutions by matching the interior solution to an exterior 
vacuum spacetime, at a junction interface, much in the spirit of [15j . 

In the analysis outlined below, we consider the graviton one loop contribution to a classical energy in a wormhole 
background [13| , where the classical energy is given by 

< - I = - 1 / ^ VgR : = - J- r - ^ 2 b '} r ] , , (5) 

7s 16ttG y s V!J 2G J ra ^l-b(r)/r r 2 g 2 (E/E P ) 1 ' 

and the background field super- hamiltonian, is integrated on a constant time hypersurface. Note that the gravi- 

ton one loop contribution to a classical energy contribution is evaluated through a variational approach with Gaussian 
trial wave functionals, and the divergences are treated with a zeta function regularization. Using a renormalization 
procedure, the finite one loop energy was considered a self-consistent source for a traversable wormhole (we refer the 
reader to [l3[ for details). 

In the following, we consider g t j — + /iy , where hij is the quantum fluctuation around the background metric 
gij . We shall also take into account the total regularized one loop energy given by 



e tt = -~y^ 91 ^ E i Ep ^ 



9 2 2 (E/E P ) 



E 2 (r) + jE 2 (r) 



(6) 
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E 2 

gl(E/E P ) h * 



where Ef (r) > 0, and Ei are the eigenvalues of 

(AIV), 

h 1 - is the traceless-transverse component of the perturbation [l6l Il7| and 

where is the Lichnerowicz operator defined by 

(A L = A/iy - 2J2i fci , fc w + R lk h) + R jk hi 

with A = — V a V Q . Rather than present all the intricate details here, we refer the reader to Ref. 
discussion on these issues. 

Using the Regge and Wheeler representation 18j, the eigenvalue equation can be reduced to 



d 2 l{l + l) 2 , ; 



d.v 



Eh 



9l (E/E P ) 



fi{x) (i = 1,2). 



(7) 

(8) 
(9) 

for a detailed 
(10) 



where we have used reduced fields of the form fa (x) = Fi (x) /r, and have defined, for simplicity, two r— dependent 
effective masses m\ (r) and m 2 (r) given by 



m 2 (r) 



ml (r) 
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7T 



3b' (r) , 3fc(r) 



2r 2 



2r 3 



(r = r (x)) 



(11) 



1 



6(r) 



fc'(r) _ 3b(r) 
2r a 2r 3 



Taking into account the WKB approximation, from Eq. (|10[) we extract two r— dependent radial wave numbers 
given by 

EL, 



1(1 + 1) 



9l (E/E P ) 



mi{r) (i = l,2). 



(12) 



It is useful to use the WKB method implemented by 't Hooft in the brick wall problem [191 ] . by counting the number 
of modes with frequency less than u>i, i = 1,2. This is given by 



g(E l ) 



Vi (l,Ei) {21 + l)dl, 



where Vi (I, Ei), i = 1,2, is the number of nodes in the mode with (I, Ei), such that (r = r (x)) 

i r +co i 

Ui(l,Ei) = - / dxJkUr,l,Ei). 

" J — OO 



(13) 



(14) 



Note that the integration with respect to x and l max is taken over those values which satisfy kf (r, I, Ei) > 0, i = 1, 2. 
With the help of Eqs. (fT3|) and (fl4| . the self-sustained traversable wormhole equation becomes 



H- 



(<>) 



1 



7T f-? 



+ °° £ 3i (E/Ep)_dHEil dE 



g 2 (E/E P ) dEi 



(15) 



The explicit evaluation of the density of states yields 



d~g{Ej) 
dE,. 



dv(l,Ei) 



1 

7T 

2 

7T 

4 
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(2/ + l)dl 
J™. (21 + 1) 



dxr 2 



o y/k*(r,l,E) dEi V ffl (S/^p) 



- < (r) dl 



dEi \g 2 (E/Ep) 



Ef 



+°° d 

dxr dEi y 9 2 2 (E/E P ) 



m l [f) 



m\ (r) 



E 



9l (E/Ep) 



m 2 (r) 



(16) 
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Finally, plugging expression (IT6l) into Eq. (fT5j) and taking into account the energy density, we obtain the following 
self-sustained equation 

1 "' (r) 2 (/ 1+ «, (17) 



2G r 2 g 2 (E/E P ) 3vr 2 

which will play an important role in the analysis below. The integrals I\ and 1 2 are defined as 

t-fBHWV^^*^-^)**, (18) 

and 



g 2 (E/E P )dE \gl{E/E P 



„ gi{E/Ep) d ( E 2 2 



g 2 (E/E P )dE \gl{E/E P 

respectively. Note that E* is the value which annihilates the argument of the root. 

In 1\ and I 2 we have included an additional 4tt factor coming from the angular integration and we have assumed 
that the effective mass does not depend on the energy E. To further proceed, we can see what happens to expression 
(fl~7|) for some specific forms of g\ (E/Ep) and g 2 (E/Ep). It is immediate to see that integrals I\ and I 2 can be solved 
when g 2 (E/Ep) = g\(E /Ep). However the classical term keeps a dependence on the function g2(E/Ep) that cannot 
be eliminated except for the simple case of g 2 (E/Ep) = 1. Therefore we will consider different models regulated 
by the Rainbow's function g\ (E/Ep), of the form given by Eq. ©, to analyse the effect on the form of the shape 
function b(r). 



III. EXAMPLES 



A. Specific case: gi (E/E P ) = expf-a-fr), ff2 (E/E P ) = 1 



Following Ref. Q, we consider the following choice for the Rainbow's functions 

g 2 [E/Ep) = 1; a > e 



gi (E/Ep) = exp(-a— 2-), 



£_ 2 

[ E\ 



Thus, the graviton contribution terms (fT5|) and pTJ|) yield the following relationships 



and 



/i = 3 



/ 2 = 3 



E I 

exp(-a-^2-)£ 2 y/£; 2 - m\ (r)dE , 



, eM~^)E 2 jE^~ml (r)dE , 



(20) 



(21) 



(22) 



respectively. Using the general results outlined in Appendix [X] for the two integrals I\ and I2, Eq. pT|) can be 
rearranged in the following way 



1 b'(r) 
2G r 2 



E% 
2tt 2 



exp 



Ki 



H exp 

a 



(23) 



where xi = ^/ mf (r) / J5p , X2 = y m-2 (r) /Sp and i^i (x) is a modified Bessel function of order 1. Note that it is 
extremely difficult to extract any useful information from this relationship, so that in the following we consider two 
regimes, namely the cis-planckian regime, where Xi <ti 1 (i = 1,2), and the trans-planckian regime, where Xj>l. 

In the cis-planckian regime, with the approximation x\ <Si 1 and x 2 -C 1, and expanding the right hand side of Eq. 
(f2"3"]). we find that the leading term is given by 



1 b'(r) 
2G r 2 



E% 
2tt 2 



\--{xl+x 2 )+0{x{ 
a 2 - a 



(24) 
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Substituting the factors Xi — \j rrif (r) / E' p (i = 1, 2) in the latter, provides the following relationship 



1 b'jr) 
2G r 2 



2ir 2 



12 



which can be rearranged to give 



Ar 2 



ar 2 Ep 



1 - 



6(r) 



2b'(r) 



aEpr 2 



b{r) \ 2b' (r) 



a 



(25) 



(26) 



where we have used the definition G = E p 2 = l 2 P . Restricting our attention to the dominant term only, we find that 



b(r) 



n 



El 
3n 2 a* 



(r 3 



(27) 



which does not represent an asymptotically flat wormholc geometry, as the condition b(r)/r — > when r —> +oo, is 
not satisfied. However, for these cases, one may in principle match these interior wormhole solutions with an exterior 
vacuum Schwarzschild spacetime. We shall not proceed with the matching analysis in this paper, but we refer the 
reader to [l5[ for further details. 

On the other hand, in the trans-planckian regime, i.e., x\ S> 1 and i2» 1, we obtain the following approximation 



_L_ V(r)_ 
2G ~7 r 



E% 



8y a a 7T 



3 -77-3 



exp (— OiXj) x\ + O 



exp (—ax%) X2 + O 



(28) 



Note that in this regime, the asymptotic expansion is dominated by the Gaussian exponential so that the quantum 
correction vanishes. Thus, the only solution is b'(r) — and consequently we have a constant shape function, namely, 
b(r) = r t . In summary, in the trans-planckian regime provided by the Rainbow's function (|20p. we verify that the 
self-sustained equation (|17|) permits asymptotically flat wormhole solutions with a constant shape function given by 
b(r) = r t . 



Specific case: gi (E/E P ) = ( 1 + ) exp(-a-f^) 



52 (E/Ep) = 1 



Another interesting choice for the Rainbow's functions is the following Q, 



g x (E/E P ) =1 + 



E 
E~p~ 



E 2 

cxp(-a— T ), 

Er 



92 (E/E P ) 



a>0,(3£ 



(29) 



For these specific functions, we once again use the general results in Appendix IA1 so that the two integrals I\ and I2, 
given by Eqs. (fT5|) - (fH)|) . take the following form 



I - E r 
h - 2 - 2^ 



a 



■ exp 



a -' 



P— [ x i,2 + ^ ) exp {-ax! 



(30) 



where once again, we have defined for notational simplicity X1.2 = ym 2 2 (r) /Ep. As in the previous example, it is 

extremely difficult to extract any useful information from these relationships, so that we consider the two regimes, 
i.e., the cis-planckian regime, where x% -C 1 (i = 1,2), and the trans-planckian regime, where x% 3> 1, respectively. 
In the cis-planckian regime, where x% 1, the self-sustained equation (|17[) takes the following form 



1 b'(r) 
2G r 2 



2air 2 



2 2y^+^f3 



2a 3 / 2 



2aV2 



(x\ + x\) + O (x 4 ) 



which leads to 



b'(r) 



E P 



4^ + 3^tt(3 2 2y/a + ^(3 ( 6 
r — 



2a 3 /2 



,1/2 



E P 



1 



b(r) 



E P 



(31) 



(32) 



Note that the term proportional to r 2 leads to a non-asymptotically flat wormhole spacetime, so that it is useful to 
do away with this term. To this effect, it is straightforward to see that for j3 = 



■ y/a/n, one arrives at 



b'{r) = 



3air 2 



(33) 
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which provides the following solution 



b(r) 



3air 2 - 12n 2 - 2 



(3a7T 2 - 2) r t 



(34) 



where we have used the condition b(r t ) — r t . It is a simple matter to verify that this solution is not asymptotically 
flat, however, as in the previous example one may match the interior solution to an exterior solution (once again, we 
refer the reader to [l5| for further details). 

On the other hand, for the trans- planckian regime, Xi 3> 1, the asymptotic series becomes 



1 b'(r) E% 



2G r 2 8VaW 



(xi + (3x1) exp (-axf) + O ( — ) + (x 2 + (3x 2 2 ) exp (-ax$) + O ( — 



(35) 



Once again, we find that the quantum corrections are exponentially suppressed, and as in the previous example leads 
to a constant shape function, b(r) = r t . 



C. Specific case: g L {E/E P ) =11+ /3#- + 7#r exp(-a 



92 (E/Ep) = 1 



Finally, consider now the following choices for the Rainbow's functions given by 



E E 



E P ' ' E 2 p 



E 1 



9l (E/E P ) = ( 1 + + 7 — ) exp(-a-^), g 2 (E/E P ) = 1. 



(36) 



Once again using the general results outlined in Appendix we find that the two integrals I\ and I 2 , given by Eqs. 
(HI])- (in]), take the following form 



h. 



/•;,', / n.r 

2^ 6XP ■ 



2 

1.2 



x\ 2 K Y 



ax\ 2 



P\fn ( n.r 



2 ) \ 1 - 2 2a 



+1*1,2 



1,2 \ , ^1,2 



1 lax 
a K A 2 



2 

1.2 



(37) 



with the definition X\,2 = \J m \ 2 ( r ) l$p- We now investigate the two limiting approximations, namely, the cis- 

planckian regime, where Xi <C 1 (i = 1, 2), and the trans-planckian regime, where xi 3> 1, respectively. 
Regarding the cis-planckian regime, the self-sustained equation (TTT)) 



1 b'(r) E% 



2G r 2 2air 2 



2 3/3^ 4 7 \ / 7 /V^F 1 



a 2 2a i / 2 a / 
2 In (axl/4) + 2 In (aa|/4) + 2 7 _e + 1 7 \ 4 

8 + 4^ ~ 167= I ^ 



0(x 6 ) 



where je is the Euler's constant. In order to simplify the analysis, consider the following imposition 



2 , 3/3 ^tt _ AS _ n 
ot? ^ 2a 5 / 2 a*— u 

j PV^ _ 1 _ n 

Q 2 2q 3 / 2 a — U 



which leads to the following solution 



(38) 



(39) 



P 



7 = —a; for a € 



(40) 



Thus, Eq. (|38J simplifies to 

1 &'(r) _ 



2G r 2 2air 2 



' in (ax 2 /4) + In (axj/A) + 2 7£ + l \ 4 

- (2^ + x 2 J 



(41) 
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Note that as the previous formula is obtained in the cis-planckian regime, the logarithmic functions in this range are 
always negative and the whole expression goes to zero from the negative side. Furthermore, as in the previous two 
examples, it is an easy matter to show that the solution is not asymptotically flat, and in principle can be matched 
to an exterior vacuum solution at a junction interface. 

Let us now fix our attention on the trans-planckian regime, where x, > 1 (i = 1,2). The integrals in the Eq. f|3T[) 
have the following asymptotic expansion 

i 12 4 ex P (-ox; 2 ) ^ (-,4 2 + ^ + (4Q - 9 ^ ^ 2 + M + 3 ^ - ^ - 15 ( 4 "" 3 7 ^ + o (x -i)) . 

' 2tt 2 PV 1 ' 2 ' a 3 / 2 y ' 1,2 4 4a 8a 32a 2 x h2 128a 3 x? 2 V 1>2; y 

Analogously to the previous cases, due to the Gaussian exponential, we find that quantum fluctuations lead to a 
vanishing contribution, independently of the choice of the parameters /3 and 7. Therefore we conclude that b'(r) = 
leading to a constant shape function, b(r) — r t . 

Collecting the results of cases described in sections Al III B|) and (jlll Cp . we conclude that the model governed 
by the Gaussian Rainbow's function g\ (E/Ep) and its polynomial variations lead only to a constant shape function 
b (r t ) — r t . In the next section, we consider a different approach to the problem of the self sustained equation by fixing 
the geometry of the wormhole. More specifically, we consider a variational approach which imposes a local analysis 
to the problem and restrict our attention to the behavior of the metric function b(r) at the wormhole throat, r t . 



IV. FIXING THE GEOMETRY 



In this section, we outline a different approach in the context of self-sustained wormholes in MDR. More specifically, 
we fix the shape function, and therefore the geometry, and restricting our analysis to the throat we find conditions 
on the specific parameter space in order to have wormhole solutions. In the semi-classical context, we emphasize that 
solutions of self-sustained wormholes were found in Refs. [l3l fl4|. by using standard regularization and renormaliza- 
tion techniques. Thus, we apply an analogous approach, however without taking into account any renormalization 
procedure, but using only the deformed spacetime at the Planck scale. 

To set the stage and for concreteness consider the specific choice for the shape function b(r) — r 2 /r [7|. Thus, from 
Eq. (fTT|) . and restricting the analysis to the throat r — r t , we find 



1 



1 



2G r 2 g 2 (E) 3tt 2 



(h + h 



with Ji and I 2 given by Eqs. (fT8|) and (fl~9|) . respectively. Maintaining the setting g^iEjEp) 
effective masses at the throat simplify to 



(43) 

1, we find that the 



m\ (r t ) = f 



3 

7T 



(44) 



Therefore I\ and I2 become 



h = 3 / 9i (E/Ep) E 2 JE 2 -^dE 

h/ri ' V r * 



(45) 



and 



ia=3 



/ g x (E/Ep) 
Jo 



E 2 JE 2 + ^dE, 



(46) 



respectively. 

Now, in order to have only one solution with variables a and r t , we demand that 



d 


1 


1 " 


d 


dr t 


~2G 


rl. 


dr t 



2 
3^ 



(h + h) 



(47) 
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which takes the following form 

3 



2 

1 = — 

7T Z 



V3/r t E P 



gi (u)u 2 Ju 2 



{r t E P y 



tcLu - 



gi (u) u 2 i u 2 + 



{rtEpY 



(48) 



E%. 



where we have set u = E/Ep and G 1 

Consider a specific form of gi(u) given by g\{u) — exp(— au 2 ) with a variable, and after integration, we find that 



_ 2 d 

7T 2 G?Q! 



^exp 



2 {r t E f 



K 



2 (r t £ P )' 



6 d 
7r 2 da 



^exp 



3a 



2 (r t S P ) z 



3a 



2 {r t E P y 



(49) 



where Ao(:r) is a modified Bessel function of order 0. Writing the explicit form of the derivatives in Eq. (|49[1 we find 

f 



f 



-f{a,x) 



2tx 2 x 2 ' 

where, for notational simplicity, we have used x = r t Ep and used the definition 
/(«.*) = ex ? (2^2) K ° (2^) " eXP (sO ^ (2^) +9eXP (-0) ^° (£) +6XP 



(50) 



0)^(^1- ( 51 ) 



3a 
2^ 



The procedure now is in principle straightforward. In order to have one and only one solution, we demand that the 
expression in the r.h.s. of Eq. (|50[) has a stationary point with respect to x which coincides with the constant value 
1. For a generic but small a, we can expand in powers of a to find 



dx 



1 



2n 2 x- 



20- 10 In (4x 2 /a) +107^ + 9 In 3 



+ O (a) , 



which has a root at 



Substituting x into Eq. (|50|) . we find 



x = r t E P = 2.973786871^- 



0.2423530631 



(52) 



(53) 



(54) 



fixing therefore a ~ 0.242. It is interesting to note that this value is very close to the value a = 1/4 used in Ref. Q 
inspired by a non-commutative analysis [201 ]. As in Ref. [14j . it is rather important to emphasize a shortcoming in 
the analysis carried in this section, mainly due to the technical difficulties encountered. Note that we have considered 
a variational approach which imposes a local analysis to the problem, namely, we have restricted our attention to the 
behavior of the metric function b(r) at the wormhole throat, r t . Despite the fact that the behavior is unknown far 
from the throat, due to the high curvature effects at or near r t , the analysis carried out in this section should extend 
to the immediate neighborhood of the wormhole throat. Nevertheless it is interesting to observe that in Ref. [l3[ the 
greatest value of the wormhole throat was fixed at r t ~ 1.16/ Ep using a regularization-renormalization scheme. From 
Eq. (|53|) , one immediately extracts rt — 1AQ/Ep which is slightly larger. 



V. SUMMARY AND DISCUSSION 



Wormhole are hypothetical tunnels that violate the null energy condition, and therefore all of the energy conditions, 
and thus it seems that a natural environment of these exotic spacetimes lies in the quantum regime, as a large number 
of quantum systems have been shown to violate the energy conditions 0, Q, such as the Casimir effect. In this 
context, it has been shown that various wormhole solutions in semi-classical gravity have been found in the literature 
[l2T - [l4l ]. In the semi-classical approach, the Einstein field equation takes the form G M „ = 8nG (T M ^) ron , where the term 
(I) 1I/ ) ren is the renormalized expectation value of the stress-energy tensor operator of the quantized field. In addition 
to this, the metric is separated into a background component, g M „ and a perturbation /i M „, i.e., g^ = <? M „ + h^. 
A key aspect is that the Einstein tensor may also be separated into a part describing the curvature due to the 
background geometry and that due to the perturbation, namely, G^ v {g a j3) — G^ v {g a fj) + AG^^, h a p) where 
AG nv{ga.p,h a p) may be considered a perturbation series in terms of /i M „. Using the semi-classical Einstein field 
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equation, in the absence of matter fields, the effective stress-energy tensor for the quantum fluctuations is given by 
8nG (T^,,) 1611 = — (AG^ u (g a fj)) Icn so that the equation governing the quantum fluctuations behaves as a backreaction 
equation. Thus, the possibility that wormhole geometries were sustained by their own quantum fluctuations, using the 
above-mentioned semi-classical approach was shown in [l3l llij . More specifically, the graviton one loop contribution 
to a classical energy in a wormhole background was taken into account, and a variational approach with Gaussian 
trial wave functionals was used. A zeta function rcgularization was involved to handle with divergences and a 
renormalization procedure was introduced and the finite one loop energy was considered as a self- consistent source 
for the traversable wormhole. 

In this work, we have considered the possibility that wormhole geometries are also sustained by their own quantum 
fluctuations, but in the context of modified dispersion relations. We considered different models regulated by the 
Rainbow's function gi (E/Ep), given by Eq. ©, to analyse the effect on the form of the shape function b(r), and 
found specific solutions for wormhole geometries in the cis-planckian regime and trans-planckian regime. In the latter 
regime, we found that the quantum correction are exponentially suppressed, thus leading to a constant shape function, 
i.e., b(r) = rt, where r* is the wormhole throat. In the cis-planckian regime, the solutions found do not represent 
asymptotically flat wormhole geometries. However, for these cases, one may in principle match these interior wormhole 
solutions with an exterior vacuum Schwarzschild spacetime. In addition to this analysis, we also fixed the geometry by 
considering the behaviour of a specific shape function b(r), through a variational approach, at the wormhole throat. 
We explored the parameter range of the Rainbow's functions and found a good agreement with previous results used in 
the literature [3| . It is important to point out that the above-mentioned variational approach presents a shortcoming 
mainly due to the technical difficulties encountered. Note that the latter variational approach considered imposes a 
local analysis to the problem, namely, we have restricted our attention to the behavior of the metric function b(r) 
at the wormhole throat, r t (this is similar to the analysis carried out in [HI)- Despite the fact that the behavior is 
unknown far from the throat, due to the high curvature effects at or near r*, the analysis carried out in this context 
should extend to the immediate neighborhood of the wormhole throat. 



Appendix A: Integrals 

In this Appendix we provide the rules to solve the integrals I\ and 1%, given by Eqs. (fT5|) - (TT^l) . with the following 
choices for the Rainbow's functions 



5i(-^) = E^^ ex P(- a w)' 92(E/E P ) = 1; tt>0,ftei (Al) 



(A2) 



rp n rpi rp2 

^/3i— exp(-a- ! 

In some particular cases such as in section IPVI the effective mass becomes negative and the integral takes the form 

f°° E 
1+ = / g 1 (—)E 2 y/E 2 +m 2 dE 
Jo &P 

Therefore we divide the integration in two classes denoted by the integrals 1+ and i_ described by 

/oo 771 
_<?i( — )E 2 y/E 2 - m 2 dE. (A3) 
E P 

For the examples discussed in section ITTT| the following relationships for i_ are extremely useful 

</-/ „ r e 2 

p 



da n 

and the following for I + 

d n I. 



_exp(-a- T )i; 2 V£ 2 ~ m 2 dE rc > 1 (A4) 
E p 



(A5) 



jn j poo t?i 

—t = (-) n E P n cxp(-a- T )S 2 " v / S 2 + w?dE n>l. 
da n Jq E p 

Taking into account gi(j^) = exp(— a-^j-) in i_, we find that changing variables E — y/x, we obtain 

I r°° x _ 

I- = — cxp(— a — ^)\fxy x — m 2 dx 

Ep fm 2 \^/3\ ( am 2 \ r/ fam 2 \ 
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where we have used the following relationship 



/>-<-V-'» P <-^ = -Lg)- roo«p(-£) W (£) *£><>. (A7) 

The same argument applies for I + to obtain 

Y roc ^ 

1+ = — / exp(— a^-)V^V x + m 2 dx 



where we have used the following relationship 



On the other hand, the integrals of the form 
d n I°_ 



can be generated by the following expression 



while integrals of the form 

poo p2 

-) n El n exp(-a-2-)£: 2 " + V£ 2 + m?dE n>2 
Jo E P 



d n I% 



da 

can be generated by the following relationship 

3 Vi 



/ dx(x + i) 1/2 xexp(-^x) = --^ + ^-^ _5/2 exp(^)(3 - 2i/z)Erfc [s/tjj] t>0,n>0. 
Jo 2 n A 



(A8) 



H {x + P) v - 1 x»- 1 e?q>{-nx)dx = ±= (^] ' V (y) exp ( ^ ) K 1/2 _ v ( ^) Rett . (A9) 

Jo V^W V 2 y V 2 y | a rg/?| < tt 



/no 
_exp^a^)E 2n+1 VE T ^^dE n>2, (A10) 



da; (x — a) 1//2 x exp (— /xx) = — ^— fi 5 ^ 2 (3 + 2/ia) exp (— /ia) a > 0, > 0, (All) 



(A12) 



(A13) 
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